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THE SPECTRUM OF A LINEARIZED 2D EULER 

OPERATOR 

Y. LATUSHKIN, Y. LI, AND M. STANISLAVOVA 



Abstract. We study the spectral properties of the linearized Eu- 
ler operator obtained by linearizing the equations of incompress- 
ible two dimensional fluid at a steady state with the vorticity that 
contains only two nonzero complex conjugate Fourier modes. We 
prove that the essential spectrum coincides with the imaginary 
axis, and give an estimate from above for the number of isolated 
nonimaginary eigenvalues. In addition, we prove that the spectral 
mapping theorem holds for the group generated by the linearized 
rS^ ' 2D Euler operator. 



i/^ ' 1. Introduction 

(N ■ 

In recent years a new interest has been drawn to understanding the 

r^ . stability spectrum of the Euler equations for the motion of inviscid fluid 

linearized about a steady state. We will not even attempt to review a 
O I vast literature on the subject, and refer the readers to the recent survey 

^ ■ [F3] and the bibliography therein, as well as to the closely related to 



c5 . this paper work in pFYj [FH| , |FVY 



^ i In this paper, we consider 2D Euler equation under periodic bound- 

> 



ary conditions. We apply Fourier transform and linearize the Euler 
equation about a steady state that contains only two nonzero complex 
conjugate Fourier modes. Two issues are addressed. 
_c^ I First, we give a full description of the spectrum of the linearization. 

Using quite different methods, some results in this direction were ob- 



tained in [pl , |L2|| . We show that the essential spectrum coincides with 
the imaginary axis, and that the discrete nonimaginary spectrum con- 
sists of finitely many points. Moreover, we give an estimate from above 
for the number of the nonimaginary isolated eigenvalues in terms of a 
transparent geometric quantity. Although the existence of nonimagi- 
nary eigenvalues is well-known, see, e.g. 0, we are not aware of any 
results in the literature that would give an estimate from above for the 
number of nonimaginary eigenvalues. The results in [ BFY | indicate 
that this estimate is sharp. 
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Second, we show that the spectral mapping theorem holds for the 
group generated by the linearized Euler operator L, that is, we prove 
that a{e*^) = e*°"^^\ t 7^ 0, for the spectrum a{-). Note that the vahdity 
of this spectral mapping property for (non-analytic) semigroups is a 
rather delicate issue. For instance, the spectral mapping property does 
not hold even for a group obtained by a first order perturbation of a 
two-dimensional wave equation, see and many other examples and 
general discussion of this phenomenon in, e.g., ||CL| , |vN|| . 

Our strategy is to use some ideas from the theory of Jacobi matrices 
and operators on spaces with indefinite metrics as well as some general 
results from the theory of strongly continuous semigroups. It can be 
summarized as follows. 

The steady state considered in this paper gives a flow on the torus 
parallel to a vector p G Z^. The linearized Euler operator L is a differ- 
ence operator acting on a space of sequences on Z^, see @. To study 
the spectrum of L we at first "slice" the grid Z^ using subsets Sq, 
q G Z^, of lines parallel to p, see (||). This gives us a way to represent 
L as a direct sum of operators Lq acting on the space of sequences on 
Z. We show that a(L) = Uq^(^q) and R{X,L) = 0gi?(A,Lq) for the 
resolvent operators. Using an appropriate rescaling and "symmetriza- 
tion", we replace the study of cr(Lq) by that of cr(-Bq), where Bq is a 
certain two-diagonal infinite matrix. The essential spectrum cress(-Bq) 
is described using Weyl's Theorem. We show that if ||q|| > ||p|| then 
-Bq is a selfadjoint operator and, therefore, Lq does not have nonimag- 
inary spectrum. If ||q|| < ||p|| then i?q is a finite rank perturbation of 
a selfadjoint operator. However, we give an appropriate choice of an 
indefinite metric that makes i?q a J-self adjoint operator on a Pontrya- 
gin space with finitely many positive squares. Standard facts about 
J-self adjoint operators give the estimate for the cardinality of the non- 
imaginary point spectrum of L in terms of the number of points in 
Z^ located inside of the open disc with the radius ||p||. Finally, the 
proof of the spectral mapping theorem is based on a general Gearhart- 
Priiss theorem, see, e.g., [vN], and an estimate for the norm of -R(A, L), 
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2. Linearization 
Consider the two dimensional Euler equations in the vorticity form: 

an 



d^ dVl du dv 
dt dx dy ' dx dy 

TT — IT is the vorticity and u 

ax ay •' 



Here f2 = curl u = ff — §^ is the vorticity and u = {u,v) is the 
velocity. We assume that u = u{x,y) and v = v{x,y) are 27r-periodic, 
x,y & [0, 27r], and have zero spatial means. If ip is the stream function, 
then u = -^,v = ^, and Q = A^p. li Q = Ekez^w ^ke^'''^"''^ 
where ci;_k = cJk, k G Z^ \ {0}, then, after a short calculation, see 
LT| , |L^] , we can rewrite (|I]) as follows: 

X^ ^(^ ^ ^' q)^k-q^q, 
qez2\{0} 

Here A{p, q), p, q G Z^, are defined by the formula 



(2) 



dt 



kGZ^ 



1 



1 



A(p,q) = A(q,p 

IPII iiq 

provided p ^ ±q, p 7^ 0, q 7^ 0, and A{p, q) : 
Here and below we denote p = (^1,^2), q 



Pi Qi 
P2 q2 



^(Q) p) = otherwise. 
(91,92), k = (A;i,fc2); 
vertical bars denote the determinant, || ■ || is the Euclidean norm. 

Fix p G Z^\{0} and F G C, and consider the steady state 00^ = 
('^k)kez2 for (0) defined as follows: 



(3) 



^^ 




where bar denotes the complex conjugate. 

If F = a + i6 and p = (pi,P2) then the vorticity Q^, corresponding 
to d;", is given by the formula 

ri^{x,y) = acos{pix + P2y) — bsm{pix + P2y)- 

The corresponding velocity field has the following components: 

u^{x,y) = (-p2asinp ■ {x,y) -p2&cosp ■ (x, y))||p||"^ 

v°{x,y) = (piasinp- {x,y) +pibcosp- (x,y))/||p||"^ 

In particular, if p = (0,j92), then the steady state (^ defines a parallel 
shear flow (see, e.g., [0, 0) with the vorticity Q^{x,y) = acosp2y — 
bsmp2y and the profile —ap2smp2y — bp2Cosp2y. For 6 = ImF = 0, 
therefore, we have the sinusoidal profile, studied in ||MS| , 0] for the case 
of a viscous shear flow and, recently, in ||BFY|, |FSV| . 
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The linearization of (H) around the steady state u;° as in 
the following operator L: 

(4) L : (cjk)kgz2 h-* ( yl(p, k - p)rcjk-p + ^(-P, k + p)ru;k+ 



gives 



kez2 



The choice of the space on which we want to consider the operator L 
is related to the choice of the space for vorticity in (p. Thus, if 1] G 
i7''(T^), r > 0, the Sobolev space, then L should be considered on the 
space l1{I?) of weighted £^(Z^)-sequences with the weight (l + Hkp)'"/^, 
k G 1?. In what follows we will consider the operator L on P'{1?), that 
is, for r = and fi G ^^(T^). 

Denote by W the shift operator W : (tUk) *-^ (^k-p)kGZ2- This is 
a unitary operator on £^(Z^). Also, consider the following operator 



p : (f^k) '— *■ (^(p, k)rco'k)5 and remark that D_p = —-Dp. 



D 

linearized Euler operator L can be represented as follows: 

L = WDp + W*D*_p = WDp - W*Dl. 

Here and below * denotes the adjoint operator, W* : (cuk) 
andD*: K) ^ (^(p, k)rcUk). 



Remark 1. The definition of A{p, q) imply that Dp 

1 



K 



Thus, the 



l^k4 



-Dp, where 



K ■■ M 



Dl 



^k 



2||pP 
1 



pi 

P2 



k2 



rtui 



2||k| 



Pi 

P2 



ki 
k2 



TiOi 



The operator D}^ is a compact operator on £^(Z^). 



O 



(iyD»k = 2-ip||-2 



rw. 



'0\* 



{ki~pi,k2-P2) 



D^*W* 



D^W 



Remark 2. Let us define the operator L° = WD^ - W*Dp*. Note that 

Pi ki - pi 
P2 k2 - P2 

implies {WD^uj)], = (D'°Vrcj)k. Therefore, {L^ 

— L", and, as a result, cr(L°) C zM. As we will see below, Weyl's 

Theorem implies that o"css(-^) = (Jess{L^)- ^ 

3. Decomposition and symmetrization 

In this section, we will first perform a "decomposition" of the oper- 
ator L. The operator L will be represented as a direct sum of certain 
operators Lq acting on the space £^(Z). Next, we will perform a "sym- 
metrization" of the operators Lq. This procedure will allow us to study, 
instead of the operators Lq, certain operators Bq that are finite rank 
perturbations of selfadjoint operators. 
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For p G Z^\{0} as above, and any q = (gi, ^2) G ^^ we denote 



(5) 



^(1 = {({ + np : n e Z}. 



Note that Sq is a subset of the hne {q + tp : t G M}, but, generally, 
Sq does not contain all points with integer coordinates that belong to 
this line. For each q G Z^ select q = q(q) such that q G Sq, ||q|| = 
inf{||q + np|| : n G Z}, and q = q + max{n : ||q|| = ||q + np||}p. The 
last condition just takes one of the two possible points q + np G Sq 
such that ||q + np|| = ||q||. Denote Q = {q(q) : q G Z^}. Clearly, 

Sq = Sg(q) and UqeQ Sii = Z2. 

Fix q G Q and let Xq = {(cuk) e ^^(Z^) : o^k = for k ^ Sq}. Note 
that Sq n Sq/ = provided q 7^ q'; therefore, {Xqjqgg is a system of 
orthogonal subspaces in £^(Z^) such that 0qgQ-^q = ^^(Z^). 

fiZ) 



Note that J : X„ 



[iU. 



q+np /n£/i 



Ci-', 



njndl 



IS an isometric 



isomorphism. Also, the operator L leaves the subspace Xq invariant. 
Thus, L = 0qgQ Lci, where Lq = VD — V*D*, and we denote 

V : (cj„) h^ (cj„_i), D = D(p, q) : (cj„) h^ (A(q + np, p)rc^„). 



Remark 3. If q = tp, t E M, then D = by the definition of A{p, q). 
Thus, in what follows we assume that q 7^ tp, that is, that Sq does 
not belong to the line {tp : t G M}. In particular, q + np 7^ and 
q + np 7^ ±p for all n G Z. O 

Our next step is to perform a symmetrization of the operator Lq. 
Fix p as above, and q G Q such that q 7^ tp, t G M. Denote 



(6) P 



2||p| 



92 P2 



in 



iPl 



|q + np\ 



neZ. 



Note that /5 7^ due to q 7^ tp. By the definition of A{p, q) we have 



A(q + np,p)r 



|q + np\ 



iPl 



Qi Pi 

q2 P2 



/3(l + 7n). 



Therefore, 

(7) Lq = {VP - V*P) diag(l + 7„)nez. 

1 _i 
Define a = i/^a//^^, \a\ = 1, and set S = aV. Note that 

Vp - V*p = -ip^S + S*)p'^ = -i\P\{S + S*). 
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Denote M^ = S + S*, observe that (M°)* = M° and remark that 



M" 












a 





a 





a 





a 












\a\ 



Let Mq = Mq diag(l + 7„)„gZ! and note that Lq = — i|/9| Mq. Thus, 

(9) a(Lq) = -2|/3|a(Mq) and L = ^ -t\(3\M^. 

qeQ 

If A ^ ^(-^)) then A ^ o'(-^q) and we have 



i?(A,L) = 0i?(A,-? 

qeQ 



|M„ 



e 

qGQ 



ri? 



M. 



The operator Mq is a multiphcative perturbation of a selfadjoint op- 
erator Mq. Instead of Mq we will consider a symmetric operator B^ 
with the same spectrum. To do this, for n G Z we define: 

(10) 6n = |1 +7n|2 if 1 + 7„ > and 5„ = i\l + 7„|2 if 1 +7„ < 0. 

Note that 6^ = 1 + 7„, and 5„ G M if and only if ||q + np|| > ||p||, and 
6n G iM if and only if ||q + np|| < ||p||. Denote 

(11) 



Eq = diag(5„)Mjdiag((5„ 



and note that 

(12) Mq = Mj dmgiSn) ■ diag((5„), B^ = diag(5„)M° diag((5„). 

Proposition 1. The nonzero elements in o'(Mq) and o'(i?q) coincide. 

Proof. This is a consequence of the following elementary fact: If A and 
B are bounded operators, then a{AB)\{0} = a{BA)\{0}. D 



" . 




" 





aS^i6o 





a6.i6o 





a6o6i 





a^o^i 






Remark 4. Assume ||q|| > ||p||. By the choice of q G Q we have 
||q + np|| > ||q|| > ||p||. Therefore, (5„ G M for all n G Z and i?q = 5*. 
Hence, (j{Bq) C M or a'(-^q) C iM. Since Lq is a bounded operator for 
each q G Q, we have that U||q||<||p|| "^(-^q) is a bounded set. Moreover, 
UqeQ o'(-^q)\'^Il^ is a bounded set. O 
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Remark 5. Assume ||q|| < ||p||. Then (5„5„+i G iM for exactly two 
values of n G Z. Each of these values corresponds to a pair of two 
consecutive points in Eq such that one of the points lies inside and 
another outside of the disc with the radius ||p||. For all other values of 
77, G Z we have SnSn+i G M. Thus, i?q is a perturbation of a selfadjoint 
operator by a rank four skew-self adjoint operator. O 

It is quite simple to describe the spectrum of the "constant coeffi- 
cients" infinite matrix M^. For q G Q and /3 = /3(q), defined in (|]), 
let 

where M° is as in (|). Note that ||M°|| < 2 and, thus, a(M°) C {\z\ < 
2}anda{Ll)c{\z\<m. 

If F : £^(Z) -^ L^(T) is the Fourier transform, then 

F((A + M°)(a;„)) = (az + A + az)F[K)](z), \z\ = 1, 

1 _i 
for each A G C, and a = i/S^/jS'^. Therefore, —A ^ cr{M^) if and 

only if the function (A + M^){-) defined by the formula (A + M^){z) = 

a'z + X + az is not equal to zero for all |2;| = 1. If this is the case, then 

(A + M^)-^ = W-^j^^^j^j^^W. Therefore, the following fact holds. 

Proposition 2. (a) a(AfO) = a^^K) = ["2' 2]/ 

We note that very interesting and deep results on the spectrum of the 
"variable coefficients" infinite matrices of the same type as Mq could 
be found in lUNII and the literature on Jacobi matrices cited therein. 



However, the results that we need for the specific rate of decay of 7„ 
to zero for the problem in hand do not seem to be available. 

4. Spectrum 

In this section, we describe the spectrum of the linearized Euler 
operator L. Let <Jp{-) denote the point spectrum. For A = a + ir, 
a 7^ 0, and /? = /?(q) defined in (||) we denote z = X/{i\ 



Theorem 3. (a) For each ci & Q we have: 

a^ULcd = cr{V(3 - V*^) = ^[-2|/?|, 2|/3|]. 

(b) a(L)=UqeQ^(^q); 

(c) acssiL) = iM. and ap{L)\i'R. = U||q||<||p|| (c^p(-^q)\'^I^) ^s o, bounded 
set with accumulation points only on zM. 
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Proof, (a). Split Lq = L^ +L™'"p, where L^ = V(3-V*f3, and L^'^p = 
Lqdiag(7„)„gz is a compact operator due to lim„^oo \ln\ = 0. Apply 
Weyl's theorem (see Lemma XIII. 4. 3 in |[RS|| ) for A = L^ and i? = Lq. 
Note that (j{A) = i[— 2|/3|, 2|/3|] has an empty interior in C, and C\a{A) 
consists of only one component. Since A — B is a compact operator, 
we have aess{L^) = o'essiLq). Using Proposition |], we have (a). 

(b). Split L into the direct sum of two operators, L* and L^, that 
correspond to "small" and "big" values of ||q||, that is, write L = 



L^ + L , where L^ 



®l|q|l<l|p|| ^"i ^^^ ^ 



01 



q||>||p|| ^q- Si^^^ ^' is 



a direct sum of finitely many operators, we have 

a{L) = { U a{L^))[ja{L'), 

l|qll<llp|| 

and we need to see only that <j{L^) C U||q||>||p|| ^(-^q)- But |/3| 
I/5(q)I ^ oo as ||q|| -^ oo. Using (a) in the theorem, we have: 

U a(Lq)D U z[-m,2 



llq||>llp|| 



l|q||>l|p|| 



Thus, it suffices to show that cr{L^) C iR. Since 

^''= ^q= -m^)\M, 

llq||>llpll llq|l>llp|l 

it suffices to prove the following claim: 



q' 



(13) 



if A ^ iM, then sup 

llq||>llp|| 



m^'^"" 



-1 



< oo. 



Indeed, assume that (|T3|) is proved. Then 



||(A-L 



b\-l\ 



m7 ^\r( 



|q||>llp|| 



«lPiqj 



(^ + Mq 



-1 



< oo 



and A ^cx{L^). 

To prove (pTS]) , fix A = a + ir, a 7^ 0. Note that for any q G Q we 



have (see (|T^)): 

2; + Mq = z + M° + Mj diag(7„) 



(14) 



OA-l /\#0 ■!• 



[z + Mj) I+{z + Ml) Mj diag(7. 
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Note that the resolvent {z + M^)~^ in this identity exists since z 
and Mq is a self-adjoint operator; moreover, 

(15) (- + M°)-^ =rr^ = fl- 

^ ^^ Im z\ \a\ 



Proposition 4. There exists a constant c(p) such that for all q G Q, 
q 7^ 0, we have: 

|/3(q)||| diag(7„)„ez|| < c(p)/||q||. 

Proof. Using (||), we have: 



diag(7„)„gz|| = \P\ sup |7„| < csup 



Qi Pi 

q2 P2 



|q + ^p 



1-2 



gi + npi pi 
q2 + np2 P2 



|q + ^p| 



csup 

nez 

csup |(q + np) ■ p"'"||q + ?2p||~^| 
nez 



where p-*- = {p2, —pi) is the Z^- vector, perpendicular to p. Using 
Cauchy-Schwartz inequality, we have that 

,-L|| 



P diag(7n)nGZ < CSUp^ ■ ^— < ^- ^ ■ : 

„6Z ||q + ^pr inf„gz||q + '^p| 

where the definition of q G Q has been used. 



llq| 



n 



For \ = a + ir as above, and c(p) from Proposition H, fix qo = qo(a) 
such that llqoll > ||p|| and if ||q|| > ||qo|| then the inequality 



(16) 



2c(p) ^ 1 
|a|||q|| ~ 2 



holds. Note that the set Q^ := {q G Q : ||q|| G [||p||, ||qo||]} is finite, 
and let Qb = {q e Q : ||q|| > ||qo||}. 

Using Proposition^, and the inequality ||Mq|| < 2, we have (thanks 
to (|1|)) that if q G Qb, then 



(^ + M0)-'Mjdiag(7 



njndl 



< — — ||diag(7^ 



<MpI<1. 

~ |a|||q|| ~ 2 



Thus, the operator I + {z + M^) ^M^ diag(7„) is invertible and, for 
qG Qb, 



(17) 



OA-i /\#0. 



/+(z + Mj) 'Mjdiag(7 



< 2. 
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For each q E Qs we remark that z ^ a{Bq), see (p!2D, and, hence, 
z ^ (j(Mq). Since Qs is a finite set, we have 

1 



sup 

qSQs 



m 



{Z + Mq) 



-1 



< oo. 



Remark 6. Let L^ = 0qgQ^ -^q- Then L^ is a bounded operator 
Hence, if A = a + ir, a 7^ 0, then ||(A-L^)~^|| =0(|r|"^) as |r| 



00. 



O 



To finish the proof of (|T|), we use (|T1-(|T5D and ([T^): 



»-ii 



sup \\\I3\-^ {z + M^) '\\< sup — 
qeQi qeQb |P| 



(z + Mj) 



o^-l 



< sup 2|/3|-i ■ ^ = A- 

This proves (|TB|) and (b) in the theorem. 
(c). Since 

0"css(-^) = ( \j 0"css(-^q))[jcrcss(-^''), 

llq||<l|p|| 

the first statement follows from (a) in the theorem. The second state- 
ment follows from Remark |^ and (a) in the theorem. D 

Since L is the sum of a skew-adjoint operator L° and a compact 
operator, see Remarks [^ and |^, the operator L generates a strongly 
continuous group. The following spectral mapping theorem holds for 
the group {e*^}i>o. Its proof is similar to the proof of Theorem 1 in 
[^JLS[. 

Theorem 5. If L is the linearized Euler operator, then 

Proof. Let V' = ®„^q^ L^- Inequality ([T8|) shows that if A = a + ir, 
a 7^ 0, then ||(A — i^^)^^|| < 2/|a|. Using Remark || we have that 

(19) \\{X-L)-^\\=0{1) as |r|^oo. 

Now the assertion follows from the resolvent estimate ( |T9|) and the 
following general Gearhart-Priiss spectral mapping theorem: On a 
Hilbert space the spectrum cr(e*'^), t 7^ 0, is the set of the points e'^* 
such that either /i„ = A + 27m/t belongs to cr{L) for some n G Z, or the 
sequence {||i?(/i„, L)||}„g2 is unbounded, see, e.g. ||CL|| . Recall, that 
the spectral mapping property always holds for the point spectrum. 
Due to (p^D we conclude that crcss(e*^) belongs to the unit circle. D 
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5. J-Theory 

In this section, we obtain an estimate from above for the number 
of non- imaginary isolated eigenvalues of the operator L. Recall that 
0"p(Lq)\2M is empty as soon as ||q|| > ||p||. Let x denote the number 
of points q G Z^ that belong to the open disk of radius ||p||, and such 
that q 7^ tp. Since for such q we have (igi)^ + (±^'2)^ < IIpII^, we 
conclude that x is even. 

Theorem 6. The number of nonimaginary eigenvalues of L (counting 
the multiplicities) does not exceed 2x. 

Proof. Since Lq = for q = tp, only those Lq for which ||q|| < ||p|| 
and q 7^ tp will contribute to the nonimaginary point spectrum of L. 
For each such q let n'^ denote the smallest and n'' the largest integer 
such that 6n € i^ for n = n'^, . . . , n'^, see (|10|) for the definition of 5„. 
Let Jq = diag(j„)„e2;, where j„ = 1 if n = rZq, . . . ,n'^ and j„ = -1 
otherwise. Note that ^q = ^q ^ = ^q and 

Jqdiag((5„) = diag(5„)Jq = (diag((5„))*. 

Using (0), we have that Jq-BqJq = -B*, that is, that the operator i?q 
is Jq-selfadjoint. Let {00,00') be the standard scalar product in £^(Z,). 
Note that the formula [co',ci;'] = {J^^uj^uj') defines an indefinite metric 
on £^(Z). Thus, (?{Ij) is a Pontryagin space with n'^ — n'^ positive 
squares. By a standard result of the theory of J-selfadjoint operators 



on Pontryagin spaces (see, e.g., [^, Cor. II. 3. 15]), we have that the 



number of nonreal eigenvalues of B^ does not exceed ^{n'' — n'). By 
(|T2|), (^), and X = Xlql'^q ~" ''^q) ^6 have the result. D 

Proposition 7. The nonimaginary eigenvalues of L are symmetric 
about the coordinate axes. 

Proof. Re-write (0) as Lq = A^diag(l + 7^), where N = V[3- V*(3 = 
—N*. Using the argument in the proof of Proposition |I], we have: 

^(^q)\{0} = ^(diag(l + 7n)iV*)\{0} = -a(diag(l + 7n)iV)\{0} 

= -a(iVdiag(l + 7n))\{0} = -a(Lq)\{0}. 



Thus, a(Lq)\{0} = (t(L*)\{0} = -a{L^)\{0}, and the nonimagi- 
nary eigenvalues of L are symmetric about iR. Next, define J on £^(Z) 
by J : {uJn) ^ ((-l)''cu„) and note that VJ = -JV and JV* = -V*J. 
Thus, 

JLqJ = J(\/diag(l + 7n)/3 - V* diag(l + 7„)^)J = -Lq, 

and cr(Lq) = a{JL(^J) = —a{Lf^). D 
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Note that the change of variables ^ = pix + p2y, r] = —p2X + piy in 
(|l|) converts Q^ to the vorticity Q^ = ^'^{C,), that is, to a parallel shear 
flow. However, this does not simplify our analysis, since the new flow 
looses 27r-periodicity, and the results from, say, cannot be applied 
directly. 

Recently, a very interesting case of a steady state whose vorticity has 
four (symmetric) nonzero Fourier modes was considered in ||FVY|| . For 



general p = (^1,^2), the steady state considered in the current paper is 
in an "intermediate position" between the parallel shear flow as in |F]] 
and the Kolmogorov flow as in [[B1^'Y[] , and the more sophisticated case 



studied in [[b'VY||. For the case in [FVT|, the continuous spectrum of 



the linearization is unstable, while in our case it is always stable (that 
is, located on the imaginary axis). Note that the vorticity Q^ for the 
case considered in [[FVYII has the following representation: 



fi"(x, y) = Re(rie'P-^^'^^ + Fae^P^'^^'^O), Fi.s G C, pi,2 G Z 
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